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The analytical probability distribution of finite systems 
obeying Bose-Einstein statistics in one, two, and three dimen- 
sions are investigated by using a canonical ensemble approach. 
Starting from the canonical partition function of the system, a 
unified approach is provided to study the probability distribu- 
tion of a condensate for various confinements and in different 
dimensions. Based on the probability distribution function, it 
is straightforward to obtain the mean ground state occupation 
number and particle number fluctuations of the condensate. 
It is found that the particle number fluctuations in a trapped 
Hose gas are strongly dependent on the type of confining po- 
tential and on the dimensionality of the system. 
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I. INTRODUCTION 

The experimental achievement of Bose-Einstein con- 
densation (BEG) in trapped dilute alkali atoms 0], spin- 
polarized hydrogen |^ and metastable helium [g| has 
stimulated considerable theoretical research on this 
unique phenomenon. The BEG has been recently real- 
ized in quasi-one and quasi-two dimensions [Q, where 
new phenomena such as quasicondensates with a fluc- 
tuating phase [q-p^ and a Tonks gas of impenetrable 
bosons p jll|jl4 rmay be observed. Among the several 
intriguing questions on the statistical properties of the 
trapped Bose gas, the probability distribution and par- 
ticle number fluctuations (^S'^Nq^ of condensate play an 
important role. Apart from the intrinsic theoretical in- 
terest, it is foreseeable that the condensate fluctuations 
will become experimentally testable in the near future 
Ip^ . On the other hand, the calculation of (iJ^A^o) is im- 
portant for investigating the phase collapse time of the 
condensate Recall that the BEG can not occur 

in one-dimensional (ID) and two-dimensional (2D) uni- 
form Bose gases at finite temperature, because in that 
case thermal fluctuations can destabilize the condensate. 
The realization of the trapped BEGs in various dimen- 
sions makes the behavior of particle number fluctuations 
a very interesting problem, especially for ID and 2D har- 
monically trapped Bose gases. 

Within a grand-canonical ensemble the particle num- 
ber fluctuations of the condensate are given by (^S'^Nq) = 



No{No-\-l), implying that SNq becomes of order N 
when the temperature approaches zero. To avoid this 
sort of unphysically large condensate fluctuations, a 
canonical (or a microcanonical) ensemble has to be used 
to investigate SNq. Because in the experiment the 
trapped atoms are cooled continuously from the sur- 
rounding, the system can be taken as being in contact 
with a heat bath but the total number of particles in the 
system is conserved. It is therefore necessary to use the 
canonical ensemble to investigate the statistical proper- 
ties of the trapped Bose gas. 

Within the canonical (or microcanonical) ensemble, 
the particle number fluctuations have been studied in the 
case of three-dimensional (3D) ideal Bose gases confined 
in a box [p^|-pO| , and in the presence of a harmonic trap 
p0|-p7[. The question of how interatomic interactions 
affect the particle number fluctuations has also been the 
object of several recent theoretical investigations [p8|-^. 
Although the phase fluctuations |^,^,|5| of a condensate 
have been studied for quasi-lD and quasi-2D Bose gases, 
to best our knowledge up to now an analytical descrip- 
tion of the probability distribution and the behavior of 
SNq for trapped ID and 2D Bose gases have not been 
given directly from the microscopic statistics of the sys- 
tem. The purpose of the present work is an attempt to 
calculate <5iVo of the trapped ideal Bose gas in various 
dimensions, using the analytical description of the prob- 
ability distribution obtained directly from the analysis of 
the canonical partition function of the system. In addi- 
tion, the probability distribution of the condensate will 
be used to calculate the mean ground state occupation 
number (A^o)- 

The paper is organized as follows. Sec. II is devoted 
to outline the canonical ensemble, which is developed to 
consider the probability distribution and particle number 
fluctuations of trapped ideal Bose gases. In Sees. III-V 
we investigate the particle number fluctuations of 3D, 2D 
and ID Bose gases trapped in a harmonic potential and 
in a box, respectively. The last section (Sec. VI) contains 
a discussion and summary of our results. 
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II. PROBABILITY DISTRIBUTION AND 
PARTICLE NUMBER FLUCTUATIONS OF THE 
CONDENSATE 

In the canonical ensemble, the canonical partition 
function of N trapped ideal bosons is given by 



Z{N)= J2 exp[-/3S„iV„e„], 



(1) 



where P = l/ksT, iV„ and £„ are occupation number 
and energy level of the state n, respectively. The basic 
starting point of our approach for calculating the proba- 
bility distribution of the condensate is to separate off the 
ground state labelled n = from the states n 7^ 0: 



N 



(2) 



No=0 



where Zq (Nt) represents the canonical partition func- 
tion of a fictitious system consisting of Nt = N — Nq 
trapped ideal non-condensed bosons: 



Zo{Nj 



E 



exp 



n^O 



(3) 



Note that the ground state n = has been separated 
off in the canonical partition function Zq (Nt) of the fic- 
titious system. In the following we will find that this 
separation plays a crucial role for the calculations of the 
probability distribution of the system. 

Assuming that Aq {N, Nq) is the free energy of the fic- 
titious system, we have 



Ao{N,No) = -kBTlnZoiNr). 



The calculation of the free energy Aq {N, Nq) is non- 
trivial because there is a requirement that the number 
of non-condensed bosons is Nt in the summation of the 
canonical partition function Zq (Nt). To proceed we de- 
fine a generating function Gq for Zq (Nt) in the following 
manner. For any complex number z, let 



Go (T, z) 



z'^^Zq{Nt). 



(5) 



The generating function can be evaluated easily. To ob- 
tain Zq (Nt) we note that by definition Zq (Nt) is the 
coefficient of in the expansion of Go (T, z) in powers 
of z. Thus one has 



Zq{Nt 



1 

2TTi 



dz 



Go (r,z) 



(6) 



where the contour of integration is a closed path in the 
complex z plane around z — Q. 



Assuming exp (z)] = Go (T, z) the saddle 
point Zq [ p6| is determined by dg {zq) / Ozq = 0. Expand- 
ing the integrand of Eq. about z = zq, we have 

Zq{Nt) = ^ /dzexp[5(zo) + o((z-zo)2)] , (7) 



where o ((z - zq)^) = i (z - ZQf f^g (zo) repre- 
sents the high order terms when expanding the integrand 
about the saddle point zo. 

Omitting the high order terms on the right hand side of 
Eq. , it is straightforward to obtain the following rela- 
tions between the free energy Aq {N, Nq) and the saddle 
point Zo of the fictitious N — Nq non-condensed bosons 



-Pj^AQ{N,No) = lnzQ. 
oNo 

In addition, the saddle point zq is determined by 



No = N-Y, 



exp [en/fcsT] Zq - 1 



(8) 



(9) 



Although the relations given by Eqs. and can 
not provide an explicit result of Zq (Nt), they are very 
useful to calculate the probability distribution of the con- 
densate. 

We should stress that the relations given by Eqs. 
and ^ are reliable although the disputable saddle- 
point method is used to investigate the free energy of 
the fictitious system. It is well known that the appli- 
cability of the standard saddle-point method for con- 
densed Bose gases has been the subject of a long de- 
bate |]l8| , ^ . In conventional approaches, the saddle- 
point method is used to discuss the canonical partition 
function Z (N) . The generating function is therefore de- 
(4) fined by G (T, z) = J2N^^^i^)- In this scheme, the 



high order terms of In Z (N) can not be omitted for the 
temperature below the onset of the BEG, because the 
factor 1/ (^ZQ^e^'^° — l) in the high order terms would be 
on the order of TV (See Eq. (8) in Ref. I^l). In our ap- 
proach, however, the saddle-point method is used only 
to discuss the canonical partition function Zq (Nt) of 
the fictitious non-condensed bosons. Because the ground 
state has been separated off in Zq {Nt), the high order 
terms of InZo {Nt) can be safely omitted. 

Using the free energy of the fictitious system, the 
canonical partition function of the system becomes 



N 



Z{N)^ Y exp[q{N,No)], 



No=0 



where 



q {N, No) = -iSNoeo - PAq {N, Nq) . 



(10) 



(11) 



Obviously, P {No\N) = exp [q {N, Nq)] /Z {N) represents 
the probability to find Nq atoms in the condensate. 
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Once q {N, Nq) is calculated from the canonical parti- 
tion function of the system, the statistical properties of 
the system can be easily obtained. However, it seems 
to be very difficult to obtain the analytical result of 
q{N,No) directly from the canonical partition function. 
To avoid this difficulty, we turn to investigate the partial 
derivative of q {N, No) with respect to iVo. Assuming 



N 



d 



dNo 

from Eq. (^l|) one obtains 



qiN,No) = aiN,No), 



(12) 



-P—Ao{N,No)=f3eo + aiN,No). (13) 

CiVQ 

Using the relations given by Eqs. (||) and (^), we get 

1 

exp [P (e„ - eo)] exp[-a [N, No)] - 1 ' 



(14) 



The most probability to find No atoms in the con- 
densate is determined by requiring gj^q {N, No) ~ 0. 
Thus the most probable value Nq is determined by set- 
ting a {N, A^o) = in the right hand side of Eq. ([l^) 



N^ = N -Y^ - 

;^oexp[/3(e„-£o)]-f 



(15) 



It is interesting to note that Nq is exactly the mean 
ground state occupation number in the frame of a grand- 
canonical ensemble. For sufficiently large N, the sum 
X]7Vo=o (0) replaced by the largest term, 

since the error omitted in doing so will be statistically 
negligible. Hence the result given by Eq. (^5|) shows the 
equivalence between the canonical and grand-canonical 
ensembles for large N. 

From the results given by Eqs. (|l|) and (|l|), we get 
the following equation for determining a (iV, No) 



no-k=y: 



1 



exp [P (e„ - eo)] - 1 



E 

n^O 



exp [13 (£n - eo)] exp[-a {N, No)] - 1 ' 



(16) 



Thus once we know the single-particle energy level of the 
system, it is straightforward to obtain a(iV, A^o). Using 
a {N, No), one can obtain the probability distribution of 
the system. 

From Eq. (p^), we get the following result for q {N, No) 

fNo 

q{N,No)= a{N,No)dNo + q{N,NP). (17) 
The partition function of the system is then 



JVn=0 



where 



G{N,No) = exp 



No 



a {N, No) dNo 



(19) 



It is obvious that G {N, No) represents the ratio 
P {No\N) /P{N^\N), i.e., the relative probability to find 
A^o atoms in the condensate. From Eq. (|l9|), the normal- 
ized probability distribution function is given by 



GniN,No)^A^e^p 



No 



a {N, No) dNo 



N^ 



(20) 



where An is a normalization constant. 

As soon as we know G {N, No), the statistical proper- 
ties of the system can be clearly described. From Eqs. 
(|l|) and (|l|) one obtains (A^o) and (S^Nq) within the 
canonical ensemble: 



(No) 



E^o=o^oexp[(?(7V,iVo)] _ ENo=oNoGiN,No) 



El=oexp[g(A^,A^o)] 



ENo=oGiN,No) 
(21) 



(S^No) = (N^) - {No) 



Eno=oNSG{N,No) 



Eno=oG(.N,No) 



Eno=oNoG{N,No) 
ENo=oGiN,No) J 



(22) 



We see that in our approach, the calculation of a {N, No) 
by using Eq. (|l6|) plays a crucial role to discuss the par- 
ticle number fluctuations of the condensate. The prob- 
ability distribution of the condensate is obtained from 
a{N,No). We can give a fairly accurate description of 
the particle number fluctuations through the calculations 
of the probability distribution, although the high order 
terms are omitted when obtaining the relations (|^) and 
(^). Our discussions on the particle number fluctuations 
are reasonable mainly due to the following reasons: 

(i) Because the ground state (n = 0) and the ex- 
cited states (n ^ 0) have been separated off when con- 
sidering the partitioin function of the system, and the 
saddle-point approximation is only used to investigate 
the canonical partition function of the fictitious non- 
condensed bosons, Zq (Nt), the high order terms omitted 
in the approximation do not give a significant correction 
to the particle number fiuctuations. 

(ii) When obtaining a{N,No) through the difference 
between A^o and A'g , the high order terms omitted in Eqs. 
(^) and (|l^) are nearly cancelled with each other. This 
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In Fig. 1 the dotted line displa ys the result of |21| , p8[ . 
Our results coincide with that of , except near the 

critical temperature. In fact, below the critical tempera- 
ture, from the probability distribution (^9|), one obtains 
the analytical result of the condensate fluctuations 



C(2) 
C(3)' 



(30) 



When obtaining Eq. (|3^) we have used Eq. (g^) de- 
rived in the Appendix. Below the critical temperature, 
Eq. @ recovers the result in |gl],|4|j2|] by noting that 
C (2) = 7r^/6. This shows the validity of the probability 
distribution given by Eq. (|2|) for discussing the statisti- 
cal properties of the system. At the critical temperature 
T3D, however, our results give 



(l-2/7r)C(2)A^ 
C(3) 



(31) 



which is much smaller than that obtained in Ref. P8[ . 
This difference is apprehensible because the analysis in 
Ref. holds in the canonical ensemble except near and 
above T^u, while our result holds also for the tempera- 
ture near T^jy. Near the critical temperature, our result 
(solid line) agrees with that of Holthaus et al. . Re- 
sults given by Eqs. (|o|) and (|3|) show that particle num- 
ber fluctuations in the 3D harmonic trap have a normal 
behavior, i.e., they are proportional to N. 

When r — > 0, the particle number fluctuations of 
the condensate should be zero. To check this, note 
that as T ^ 0, G3D {N,No) = A^^ if A^o = N, while 
GsD {N, No) when Nq 7^ A^. Additionally, N 
in the case of T — > 0. Thus, we obtain {Nq) N and 
(S^Nq) when T ^ 0. The correct description of 5No 
near zero termpcrature and critical temperature again 
shows the validity of our method. 

It is well known that confinement reduces the fluctu- 
ation effect. Thus the fluctuations for bosons confined 
in a box should have a stronger dependence on A^, in 
comparison with the case confined in the harmonic trap. 
To show this we consider a 3D ideal Bose gas confined 
in a cubic box. The single-particle energy level in this 



case is given by 



/2rnL^, where 



L is the length of the box. From Eqs. ( |16| ) and (19), the 
distribution function of the condensate is given by 



G3B {N,No) ^ A^geyi-p 



127rV"2 



1^0 - A^o^l' 



(32) 



where V ^ L^, A^b and A = ^^2^^ l31r/ m are normal- 
ization constant and thermal wavelength, respectively. 
From Eq. (|3^) we obtain the particle number fluctua- 
tions below the critical temperature: 



(33) 



where the coefficient A = 2/77'' x E„^ol/n'' = 0.105. We 
see that, different from the harmonic trap, the particle 
number fluctuations of the 3D Bose gas in a box exhibit 
an anomalous behavior, i.e., it is proportional to y^/^ 
(or Af4/3)_ 



IV. TWO-DIMENSIONAL BOSE GASES 

In this section we turn to discuss trapped 2D Bose 
gases. The single-particle energy level of a 2D harmoni- 
cally confined ideal Bose gas takes the form 



1 



ny + -] hWy. 



(34) 



From Eq. ( |l4| ) one has 



E 



A^o = A- 
1 



exp [P {uxhux + Uyhuy)] exp[-a (A^, A^o)] - 1 



(35) 



Then we obtain Aq through a simple integration over n 
0: 



A^o A- 



huj' 



2D 



92 (a) 



[({2)]^/^ ksT In N 



2D 



(36) 



where W2D = i^x'^y)^^'^- When < a << 1, there is a 
good approximation for g2 (a): 

52(a) «C (2) + a (1- In a). (37) 

In the case of a < and \a\ << 1, 52 (a) is 

52(a)«C(2)-|a|(l-ln|a|). (38) 



ble value A'q as: 



By setting a = in Eq. (36), we obtain the most proba- 



/ ksT 



N^^N 

\ Tllu 



2D 



C(2) 



[C(2)]i/2fc^rin7V 



2D 



(39) 



Using the approximations given by Eqs. ( |37| ) and (|3q), 
we arrive at the following result for a {N, Nq): 



a{N,NQ) 



C{2)C{3){Nq-NS) 
2Nt^ In {Nf^/C (2)) ' 



(40) 



/ \ 1/2 

where T2D — (^^j ^t|^ critical temperature 

of the system in the thermodynamic limit. Accordingly, 
the probability distribution of the system reads 



G2Z)(A,Ao) = A2^exp 



C(2)C(3) (Nq-N^)' 
mt^ In {Nt^/C {2)) 



(41) 
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APPENDIX 

The probability distribution of a harmonically trapped 
Bose gas is a Gaussian function for various dimensions. 
The ratio between P{No\N) and P{Nl^\N) is given by 



G{No,NP)=exp -a{No-NP) 



(57) 



where a is a constant and is related to the confinement, 
dimensionality and the total number of particles of the 
system. The sum in Eq. (^2]) can be replaced by an 
integral 



N 



NoG {N, No) dNo 



No=0 



{x + Nq) exp[—ax^]dx. 

(62) 



Using the formulas (|60|), ( |6l| ) and (|2|), it is straight- 
forward to obtain the analytical result of (^S^Nq') for the 
temperature below the critical temperature 



1 

2a' 



(63) 



We can see from Eq. (|63|) that the behavior of the fluc- 
tuations is determined by the factor a. 

At the critical temperature in the thermodynamic 
limit, we have A^q — 0. The probability distribution of 
the system is then 



G{No,NP) = exp [-aN^] 



(64) 



After a simple calculation, the analytical result of 
at the critical temperature is given by 



(<5'Wo> 



TT 2a 



(65) 



/^„=0^oG(iV,iVo)ciiVo 



lN,^oG{N,No)dNo 



/j(I=o^oG(iV,7Vo)diVo 



J^^^^G{N,No)dNo 



(58) 



Using the integral transformation x — Nq — Nq, one 
obtains 

NlG {N, No) dNo = / {x + N^)^ eyip[~ax'^]dx. 

No=Q J-Ng 

(59) 

For the temperature below the critical temperature, 
A^Q >> 1. Hence the integral in Eq. ( ^9|) can be es- 
timated to be 



N 



No=0 

Similarly, we have 



N^G {N, No) dNo ~ {x + N^)^ exp[-ax^]dx. 



(60) 



/ G{N,No)dNof=i / cxp[-ax^]dx, (61) 

Jno=0 J-ao 



and 
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